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ABSTRACT 

Spectrophotometry in the AA 3400-7400 range is presented for 13 areas of the bright- 
est H II region in the SMC: NGC 346. The observations were obtained at CTIO with 
the 4-m telescope. Based on these observations its chemical composition is derived. 
The helium and oxygen abundances by mass are given by: Y(SMC)= 0.2405 ± 0.0018 
and 0(SMC)= 0.00171 ± 0.00025. From models and observations of irregular and blue 
compact galaxies it is found that AY/ AO = 3.5 ± 0.9 and consequently that the pri- 
mordial helium abundance by mass is given by: Y p = 0.2345 ± 0.0026(lcr). This result 
is compared with values derived from Big Bang nucleosynthesis, and with other deter- 
minations of Y p . 

Subject headings: galaxies: abundances — galaxies: individual (SMC) — galaxies: ISM — 
H II regions — ISM: abundances 

1. Introduction 



The determination of Y p based on the Small Magellanic Cloud can have at least four significant 
advantages and one disadvantage with respect to those based on distant H n region complexes: a) 
no underlying absorption correction for the helium lines is needed because the ionizing stars can 
be excluded from the observing slit, b) the determination of the helium ionization correction factor 
can be estimated by observing different lines of sight of a given H n region, c) the accuracy of 
the determination can be estimated by comparing the results derived from different points in a 
given H II region, d) the electron temperature is generally smaller than those of metal poorer H II 
regions reducing the effect of collisional excitation from the metastable 2 3 S level of He I, and e) the 
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disadvantage is that the correction due to the chemical evolution of the SMC is in general larger 
than for the other systems. 

The determination of the pregalactic, or primordial, helium abundance by mass Y p is paramount 
for the study of cosmology, the physics of elementary particles, and the chemical evolution of galax- 
ies (e. g. Fields & Olive 1998; Izotov et al. 1999; Peimbert & Torres-Peimbert 1999, and references 
therein). In this paper we present a new determination of Y p based on observations of the SMC. 
This determination is compared with those carried out earlier based on extremely metal poor 
extragalactic H n regions. 

2. Observations 

Long slit spectra were obtained at CTIO during two observing runs, in August and September 
1990, with the 4-m telescope equipped with the R-C Spectrograph and a coated GEC CCD detector. 
Using three different gratings (at first order) the spectral ranges between AA 3440-5110, 4220-7360, 
and 5800-7370 were covered. The slit, oriented E-W, was 4.7 arcmin long and 1.6 arcsec wide; 
the scale along the dispersion axis was 0.73" /pixel. The resolution was 7 A for the blue and red 
wavelength ranges and 14 A for the intermediate one. The slit was placed at two different positions 
of the nebula, five extraction windows were defined in one slit position (regions 1-5) and eight 
in the other (regions 11-18), four of the extraction windows included the brightest ionizing stars 
(m~14), while the other 9 positions avoided stars brighter than m = 17 to minimize the stellar 
contamination of the nebular spectra. Table 1 presents the positions and sizes of the extraction 
windows. 

In Table 1 we also present the observed H/3 fluxes, F(Hf3), before correction for extinction. 
Massey, Parker, & Garmany (1989) present excellent pictures of the stellar cluster of NGC 346 
where the ionizing stars included in Table 1 are indicated, and the observed regions can be located. 
Ye, Turtle, & Kennicutt (1991) also present excellent pictures of NGC 346 where the filamentary 
structure of the nebula can be appreciated. Also in Table 1 we define regions A and B. Region 
A is the sum of regions 2, 3, 5, 12, 13, 17, and 18. Region A was defined to minimize the errors 
of the emission lines and the effect of stellar underlying absorption; the main results of this paper 
will be based on it. Regions 1 and 11 were not included in region A because they are fainter than 
the others (see Table 1) and the equivalent widths of their emission lines are also smaller (see the 
last paragraph of this section), both effects could increase the systematic errors. Region B is the 
sum of all observed regions and will be used to show the effect of the underlying absorption on the 
measured line intensities. 

To flux calibrate the spectra five or six spectrophotometric standards, from Stone & Baldwin 
(1983), were observed each night with the slit widened to 6.4". A He-Ne-Ar lamp was used to 
perform the wavelength calibration. Dome-flats and sky-flats were obtained to flatten the red 
frames while a quartz lamp flat and a sky flat were used for the blue range frames. Several bias 
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frames were taken each night. Data reduction was performed using the IRAF reduction package. 
Sky subtraction was made from observations taken one degree away from the nebula before and 
after each nebular observation. Short exposure frames at each position were used to remove cosmic 
rays and to measure fluxes of strong lines saturated in the longer exposures. 

In Figures 1, 2 and 3 we show spectra of region A at different wavelength ranges. In Figure 4 
we compare regions A and B near the Balmer limit to show the effect of the underlying absorption 
on the Balmer lines and on A 4026 of He I. 

In Tables 2, 3, and 4 we present the intrinsic line intensities, /(A), given by 



where F(X) is the observed line flux corrected for atmospheric extinction and C(H/3) is the log- 
arithmic reddening correction at H/?, and /(A) is the reddening function. For /(A) we adopted 
the normal extinction law (Whitford 1958). C(R/3) was obtained by fitting the observed Balmer 
decrement with that computed by Brocklehurst (1971) for T e = 12 000 K, and N e = 100 cm~ 3 ; the 
Balmer decrement is almost insensitive to the expected variations in T e and N e over the observed 
volumes. Table 2 presents the line intensities for regions A and B that were obtained after adding 
the spectra of all their components, while Tables 3 and 4 present the line intensities for regions 1-5 
and 11-18 respectively. 

The C(H/3) value for region A amounts to 0.15 ± 0.01, a value in good agreement with the 
values derived from the stellar data of the cluster by Niemela, Marraco, & Cabanne (1986) and 
Massey et al. (1989) that amount to 0.20 and 0.18 respectively. We have adopted the C(H/3) value of 
region A for all the observed regions. The differences between the adopted /(A) and the reddening 
function in the visual region in the direction of NGC 346 are small, moreover their effect on the 
determination of the line intensity ratios is negligible because with the adopted C(H/3) value we 
have recovered the theoretical Balmer decrement of the four brightest lines that are not expected 
to be affected by underlying absorption. 

The SMC reddening law in the UV is very different to the normal Galactic one, but fortunately 
in the visual both laws are similar (Bouchet et al. 1985); according to Bouchet et al. (1985) R = 
Ay /E(B — V) amounts to 2.7 ± 0.2 for the SMC. From the use of various reddening laws with R 
in the 3.0 to 3.2 range (Whitford 1958; Nandy et al. 1975; Seaton 1979) we estimate that the error 
introduced by the adopted reddening law is about 0.002 dex for all line ratios and smaller than 
0.001 dex for ratios of lines closer than 500 A. 

In Tables 2, 3, and 4 one-cr errors are included, the errors were estimated by comparing the 
results derived from the two different observing seasons. The total number of photons received 
by each of the He I lines of region A is in the 3 x 10 4 to 6 x 10 5 range, therefore the errors 
presented in Table 2 are about a factor of five larger than those given by photon statistics. The 
errors presented for region A are larger because they include other sources of error present in the 
reduction procedure. 
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The line intensities in Tables 2, 3, and 4 were not corrected for underlying stellar absorption 
with the exception of the Ha lines of regions B, 4, 14, 15, and 16, where we adopted an underlying 
absorption of 2 A; for these regions we adopted the theoretical I (Ra) / 1 (R/3) value to normalize 
all the lines to Ha, notice that H/3 has an intensity smaller than unity because it has not been 
corrected for underlying absorption. 

The A 3727 line intensity was corrected for the contribution due to H13 and H14 and to A 3724 
of [S III], these contributions were estimated from the other Balmer lines and from the A 6312 of 
[S III]; 7(4711) of [Ar iv] was obtained after subtracting the expected contribution of A 4713 of He I 
based on the work by Smits (1996); 7(3889) of He I was obtained after subtracting the expected 
contribution of H8 based on the work by Brocklehurst (1971). 

In Table 5 we present the H I and He I equivalent widths for all the observed regions. The 
equivalent widths of regions B, 4, 14, 15, and 16 are strongly affected by the stellar underlying 
absorption. No correction due to underlying absorption has been made for any line in this table. 
After correcting region A for extinction, based on the four brightest Balmer lines, it is found 
that the weaker Balmer lines (H9 to H12) are not affected by stellar underlying absorption (see 
Figure 4), therefore the He lines are not expected to be affected by underlying absorption. From 
similar arguments we expect each of the regions included in region A (2, 3, 5, 12, 13, 17, and 18) 
to be unaffected by underlying absorption. 



3. Temperatures and Densities 

We derived the temperatures and densities presented in Tables 6 and 7 based on the program 
of Shaw & Dufour (1995) for forbidden lines. T(0 in), T(0 n), and N e (S n) values were derived 
from the 4363/5007, 3727/7325 and 6716/6731 ratios respectively; to derive T(0 n) we considered 
the contribution by recombination to the AA 7320,7330 line intensities (Liu et al. 2000). We will 
define another temperature given by 

T(0 „ + O m) = JV(0+)n0 ' l) JV + ^ ) (0++)n0m) . (2) 

where N(0 + ) and N(0 ++ ) are obtained in Section 4 (see Table 7). 

To derive the root mean square density, iV e (rms), we adopted the following equation 



,'3<P\ ( /(Ha) \T A(He+) iV(He ++ )" 1 
Norms' - ' ] 1 ] ' 1 



'e V ••»"> I 3 



a(Ra)hu(Ra) 



1 H + 2 

JV(H+) JV(H+) 



(3) 



where a(Ha) is the effective recombination coefficient (e. g. Brocklehurst 1971), d is the distance to 
the SMC (64 kpc Reid 1999, and references therein), and r is the radius of the adopted homogeneous 
sphere (150 arcsec). From the He + /H + and He ++ /H + values for region A presented in the next 
section, the a(Ha) value for T = 11950 K, and the I(Ha) measured by Kennicut & Hodge (1986) 
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we find that iV e (rms) = 14cm -3 (see Table 7). It can be shown that, in the presence of density 
fluctuations (which is always the case), JV e (rms) provides us with a lower limit for the local density, 
A^ e (local). To derive the abundance ratios we need to use N e (local) values; in the presence of 
density fluctuations A" 2 (rms) = e N 2 (local), where e is the filling factor. From our iV e (rms) value 
and the iV e (He 

ll)sc present in Table 7 (defined in the last paragraph of this section) it follows 
that 0.01 is a representative value for e. 

From the ratio of the Balmer continuum flux to a Balmer line flux it is possible to derive 
the temperature T e (Bac), from the Balmer line emissivities (Storey & Hummer 1995; Hummer Sz 
Storey 1987; Brocklehurst 1971) and the continuum emissivities for the He I and H I continuum 
(Brown & Mathews 1970) we find the T e (Bac) value presented in Table 7. 

By combining T(Bac), T(0 II), and T(0 in), together with N(0 + ) and N(0 ++ ) (see Section 4), 
and assuming that t 2 (0 n) = t 2 (0 III), it is possible to determine the mean square temperature 
variation i 2 (H n) = t 2 , the average temperature To(H n) = To, and t 2 (0 III) over the observed 
volume since (Peimbert 1967) 

. _ fT e N e N(X+)dV 
To(X } " jN e N(Xi+)dV ' (4) 

,2 (xi+) = J(r e -To(X i +)) 2 jV e jV(X i +)^ 
1 ' T (X i +) 2 /Ar e Ar(X i +)dF ' { ' 

T e (Bac) = T (l - 1.67t 2 ), (6) 
T(0 m) = T (O m) 



( 90800 _ \ t 2 (0 m) 
1+ U(0 hi)" " 



2 



T(0 n) = T (O n) 



1 + 



/ 97300 o ^ t 2 (0 n) 
\T (O ii) " 



2 



(7) 
(8) 



jV(O+)T (O ii)+iV(0++)To(0 m) 
To = N(0) • (9) 

In region A T(0 n) < T(0 ill), therefore T (O n) / T (O in), and from equation (5) it follows 
that t 2 > t 2 (0 in). In Table 7, we present the t 2 (Bac, O n+O in) value for region A. 

In the low density and low optical depth limit it follows that the emissivities of the helium 
and hydrogen lines are proportional to powers of the temperature and consequently that T e (He n) 
is given by (Peimbert 1967) 

T e (Heii) = T e (Heii,Hn) = T [l + ((a) + (3 - l)t 2 /2] 

= T (l - 1.43t 2 ), (10) 

where (a) is the average value of the power of the temperature for the helium lines and (5 for H/3. 
The a powers in the low density limit for the AA 3889, 4026, 4388, 4471, 4922, 5876, 6678, 7065, 
and 7281 lines are -0.72, -0.98, -1.00, -1.02, -1.04, -1.12, -1.14, -0.55, and -0.60, respectively (Smits 
1996), and (3 = —0.89 (e. g. Brocklehurst 1971). In the low density limit and weighted according to 
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the observational errors we obtain (a) = —0.96; for iV e (He n) = 143 cm -3 we obtain (a) = —0.89. 
From i 2 (Bac, O n+O in) and To(Bac, O n+O in) derived from equations (6 - 9) together with 
equation (10) we obtain T e (He u)boc = 11890, which is a representative temperature for the He I 
lines. 

Based on nine emission line ratios in the next section we derive A r (He + )/iV(H + ), T e (He n) 
and iV e (He II) self-consistently, hereinafter iV e (He ll)sc and T e (He n) sc (which are presented in 
Table 7). In Table 7 we also include i 2 (He II, O n+O in) derived from equations (7), (8), (9), and 
(10); and i 2 (Bac, He II, O n+O in), an average of the two t 2 determinations; T e ((He n)) which is 
the average of T e (He n) Bac and T e (He i\)sc ; and N e ( (He n) ) which is the density that corresponds 
toT e ((He II)). 

From this discussion and the values in Table 7 we conclude that t 2 = t 2 (H n) = i 2 (He n) = 
0.022 ± 0.008 and that to a very good approximation i 2 (0 m) = t 2 (H n). 

4. Ionic Chemical Abundances 

To determine the abundances from collisionally excited lines many authors adopt a two tem- 
perature scheme, with T(0 III) for the high degree of ionization zones and T(0 II) or T(N n) for 
the low degree of ionization zones. Therefore to compare with the abundances determined from 
collisionally excited lines by other authors we will assume also that that the temperature within 
the O II and O in zones is constant, that is t 2 (0 n)=t 2 (0 m)=0.000. Two points should be stated 
here: (1) under the assumption that t 2 (0 n)=t 2 (0 m)=0.000 and since T(0 II)^ T(0 III), from 
equation (5), it follows that t 2 ^ 0.000; from the T(0 in) and T(0 n) values for region A (see 
Table 6) and equations (5) and (9), we find that t 2 = 0.0013; (2) the t 2 = 0.0013 abundances are 
a lower limit to the real abundances since the assumption that t 2 (0 n)=t 2 (0 in)=0.000 implies 
constant temperature within the O n and O in zones, which is not the case (see section 3). The 
total abundances for i 2 (H n)=0.022, our preferred t 2 value, will be discussed in section 5. 

Therefore we have determined the ionic abundances of the heavy elements for all the regions 
observed using t 2 = 0.0013, the abundances are presented in Table 8. The abundances were 
computed with the program presented by Shaw & Dufour (1995). To determine the ++ , Ne ++ , 
S ++ , Ar ++ , and Ar 3+ abundances we used T(0 III), while for N + , + , and S + we used T(0 II). 
For T(0 in) we adopted the value derived for each region (see Table 6). On the other hand, it 
can be seen in Table 6 that the T(0 n) temperatures are systematically higher for those regions 
that include bright stars and have relatively smaller EW(H.(3) values (regions 4, 14, 15, 16, and B) 
than for the regions that have larger EW(H/3) values; however, this effect is not real and has to do 
with the difficulty of establishing a proper continuum baseline owing to the presence of underlying 
Balmer lines in absorption in the 3700-3750 A region, where the A3727 lines originate. This effect, 
together with large T(0 n) errors presented in Table 6 , which are due to other causes, led us to 
assume for all regions that T(0 n)= 0.9036 T(0 in), the value determined for region A. 
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To estimate the B./3 emissivity we adopted 



/(H/3) = i (up n 0+) + /(H/3 n o ++ ), 

/(H/3nO+) _ T(Oii)-°- 89 JV(0+) 
/(H/3nO++) ~ T(0 in)" - 89 iV(0++) ' 



(11) 



(12) 



where /(H/3nO + ) and I(H/?nO ++ ) are the intensities of H/3 in the + and ++ regions respectively. 
We consider our procedure adequate since the total O abundances derived for all the regions 
show a small dispersion despite their different ionization degree (compare for example the total O 
abundances of regions 16 and 18 that are practically the same). 

To obtain He + /H + values we need a set of effective recombination coefficients for the He and 
H lines, the contribution due to collisional excitation to the helium line intensities, and an estimate 
of the optical depth effects for the helium lines. The recombination coefficients that we used were 
those by Storey & Hummer (1995) for H, and Smits (1996) for He. The collisional contribution 
was estimated from Kingdon & Ferland (1995) and Benjamin, Skillman, &; Smits (1999). From 
the intensity of A3614 and the computations of Robbins & Bernat (1973) it was found that the 
He I singlet lines were produced under case B. The optical depth effects in the triplet lines were 
estimated from the computations by Robbins (1968). 

To derive the He + /H + value of region A, in addition to the Balmer lines, we made use of nine 
He I lines, AA 3889, 4026, 4387, 4471, 4922, 5876, 6678, 7065, and 7281 to determine iV e (He n) and 
T e (He n) self-consistently. Each of the 9 He i/H I line ratios depends on T e , N e , N(Re + )/N(H+), 
and the optical depth of 3889 (T3889), and each dependence is unique. Therefore we have a system of 
9 equations and 4 unknowns. We decided to obtain the best value for the 4 unknowns by minimizing 
X 2 - The x 2 value is given by: 



where cr(Aj) is the absolute error in the measurement that can be obtained from Table 2, and 
A^(He + (Aj, T e , A^e, T3889))/-/V"(H + ) is the abundance derived from each line for those parameters. 
The best T3889 value is slightly negative which is unphysical, moreover from CLOUDY models 
(Peimbert, Peimbert, & Luridiana 2000) it is found that T3889 is close to zero, therefore we decided 
to adopt T3889 = 0.0 and to use equation (13) to derive T e , N e , and N(Re + )/N(R + ) self-consistently. 
For a system with 9 independent determinations and 3 unknowns we have 6 degrees of freedom and 
we expect the minimum x 2 to be in the range 1.64 < x m in < 12.59 at the 90% confidence level. 
The value of Xmin = 6-53 found in Table 9 is in excellent agreement with this range. In Table 10 
we present He + /H + values for different temperatures and densities; the temperatures were selected 
to include T(0 III), T(Bac), T(He u)sc, and two representative temperatures; the densities were 
selected to include the minimum x 2 at each one of the five temperatures. The temperature with 
the minimum x 2 is the self-consistent T(He 11) and amounts to 11950 ± 560 K; this temperature is 
in excellent agreement with the temperature derived from the Balmer continuum that amounts to 




(13) 
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11800 ± 500 K, alternatively T(0 III) amounts to 13070 ± 100 K. Notice that the \ 2 test requires 
a higher density for a lower temperature, increasing the dependence on the temperature of the 
He + /H + ratio. As mentioned above, the values in Table 9 correspond to the case where T3889 
equals zero, for higher values of T3889 the \ 2 values increase. In Table 9 we present He + /H + values 
for a set of temperatures and densities. 

From Table 9 we obtain that T e (He 11) S c = H950 K and iV e (He 11) S c = 143 cm' 3 , which 
correspond to He+/H+ = 0.0793. 

From the surface of T e , N e and He + /H + values defined by the condition \ 2 = Xmin + 1 the 
one-cx errors presented in Table 7 were computed. 

In Table 10 we present the He + /H + values for all the observed regions without underlying 
stellar absorption. We have adopted the iV e (He ll)sc and T e (He ll)sc values of region A to deter- 
mine the He + /H + values for all the other regions in Table 10 because their observational errors are 
higher than for region A and do not permit to obtain A e (He 11) and T e (He 11) self-consistently for 
each region. 



5. Total Abundances 

In Table 11 we present the total abundances of region A for t 2 (0 n)=t 2 (0 m)=0.000 and 
t 2 = 0.0013. The gaseous abundances were obtained from the following equations (Peimbert & 
Costero 1969) 

N(Q) N(0+) + N(0++) 

7V(H) N(R+) ' 1 ' 

N(N) / A(Q+) + A(Q++) \ iV(N+) 

JV(H) V N(0+) ;at(h+)' 1 ' 

TV(Ne) {N(0 + ) + N(0 ++ )\ N(Ne ++ ) , N 

- ' ' (16) 



A(H) V N(0++) J A(H+) ' 

The abundances derived for all the other regions based on these equations agree within the 
errors with those derived from region A indicating the reliability of the equations. To obtain the 
total O abundances we assumed a correction of 0.04 dex due to the fraction of O tied up in dust 
grains, this fraction was estimated from the Si/O values derived for the Orion nebula (Esteban 
et al. 1998) and for H 11 regions in irregular galaxies (Garnett et al. 1995). 

To obtain the total S abundance we adopted the following equation 

A(S) _ rrF(Q , N(S+)+N(S++) _ 

from Table 8 it can be seen that the higher the O ionization degree the lower the [iV(S + ) + 
iV(S ++ )]/A r (H + ) ratio indicating, as expected, the increase of S 3+ with the O ionization degree. 
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Therefore to obtain the S abundance we decided to take regions 11 and 12 as representative of 
NGC 346 since they are expected to have the smallest ICF(S) values. Therefore from the ICF(S) 
values computed by Garnett (1989), the data for regions 11 and 12 in Table 8, and equation (17) 
we derived the value presented in Table 11. 

To obtain the total Ar abundance we adopted the following equation 

im ~ icf{At) tv(h+) ' (18) 

where ICF(At) includes the Ar + /H + contribution and according to Liu et al. (2000) can be ap- 
proximated by 

ICF(Ar) = (l- N A i°^ ) • (19) 



N(0) 

Based on the T e values presented in Table 7 we also include in Table 11 the total abundances 
for t 2 = 0.022, our preferred t 2 value, following the procedure outlined by Peimbert & Costero 
(1969). The relevant equations are: 



N{X i+ ) T{Rl3{X i+ })-°- 89 T(X,X t+ y 



-0.5 



iV(H+) T(A/A',X i +)- 1 - 39 

/ AE AE 



x exp 



VfcT(A,X*+) fcT(A/A',X i +) 
N{X i+ ) 



_N{R+) 



(20) 

A/A' 



where AE is the energy difference between the ground and the excited levels, A is the wavelength 
of the nebular line (5007 A for O in), A' is the wavelength of the auroral line (4363 for O in), 
T(H/3{X* + }) is the T(H/3) value in the region where the X i+ ion is present and is given by 

T(H/3{X i+ }) = T (X i+ ) (1 - 0M5t 2 (X i+ )) ; (21) 

T(A,X i+ ) is the representative temperature of the nebular line which can be obtained from 

T(X,X i+ ) = T (X i+ ) 

f \ (AE/kT (X^)) 2 - 3Aff/fcr (X»+) + 3/4 ] t 2 (X>+) \ 

\ + [ AE/kT (Xi+)-l/2 \ 2 )' Kll) 

and [N(X i+ )/N(R + )] x/xl is the abundance derived from T(A/A') (T(0 III) for O III). 

It is possible to derive the abundances for other t 2 values by interpolating or extrapolating the 
abundances presented in Table 11. 

In Table 11 we compare the abundances of NGC 346 with an average of the H n regions in the 
SMC derived by Dufour (1984) and those of the Sun and M17. Notice that the abundances derived 
by Dufour should be compared with those for t 2 = 0.0013, and that no correction for the fraction 
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of O tied up in dust grains was included by Dufour. Alternatively the M17 values are for t 2 = 0.037 
and 0.08 dex has been added to the O abundance to correct for the fraction of this element tied 
up in dust grains, consequently the M17 values can be compared directly with the NGC 346 values 
for t 2 = 0.022. 

The total He/H value is given by: 

iV(He) N (He ) + iV(He + ) + 7V(He++) 



7V(H) iV(H°) + iV(H+) 

The He ++ /H + ratio can be obtained directly from the 4686/H/3 intensity ratio. In objects 
of low degree of ionization the presence of neutral helium inside the H II region is important and 
ICF(He) becomes larger than 1. The /CF(He) can be estimated by observing a given nebula at 
different lines of sight since He is expected to be located in the outer regions. Another way to deal 
with this problem is to observe H n regions of high degree of ionization where the He amount is 
expected to be negligible. 

Vilchez & Pagel (1988, see also Pagel et al. 1992) defined a radiation softness parameter given 

JV(Q+)JV(S++) , 
^ N{$+)N{0++Y { ' 

for large values of £ the amount of neutral helium is significant, while for low values of £ it is 
negligible, where the critical value is around 8. In Table 10 we present the ( values for all the 
observed regions, the £ values indicate that the amount of He inside the H + region is negligible. 

On the other hand, for ionization bounded objects of very high degree of ionization the amount 
of H° inside the He + Stromgren sphere becomes significant and the ICF(He) can become smaller 
than 1. This possibility was firstly mentioned by Shields (1974) and studied extensively by Armour 
et al. (1999), Viegas, Gruenwald, & Steigman (2000), and Ballantyne, Ferland, & Martin (2000). 

To study this possibility we have estimated the ICF(He) from three different methods. (1) 
We have divided the regions in Table 10 in three groups: H (regions 5, 17, and 18), I (regions 1, 2, 
3, and 13) and L (regions 11 and 12), where H, I, and L stand for high, intermediate and low degree 
of ionization regions; within the errors the three groups yield the same He/H ratio, indicating that 
the JC-F(He) is very close to unity. (2) Ballantyne et al. (2000) have defined the following cutoff 

([O m] A5007/H/3) cutoff = (1.139 ± 0.306) + (2.5 ± 0.4)O/H x 10 4 . (25) 

For observed values higher than the cutoff ICF(He) is very close to unity; from Table 11, and 
without considering the fraction of O embedded in dust grains, we have that log N(0)/N(R)+12 = 
8.11 and, consequently, a cutoff of 4.36 ± 0.60. Alternatively the observed value is 5.43 ± 0.03. (3) 
Ballantyne et al. (2000) also find that for [O m]A5007/[O i]A6300 > 300, the JCF(He) becomes 
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very close to unity; from Table 2 it is found that 7(5007) //(6300) > 600. (It should be noted that 
the [O i] lines present in Figure 2 are blends of telluric and nebular lines even though we subtracted 
the sky contribution. Unfortunatelly the [O i] line intensities in the sky varied in a scale of minutes 
leaving a telluric remnant present, which becomes apparent because the centroids of the [O i] lines 
are blue shifted by about 3 A from the centroid of the other lines.) From these three methods we 
conclude that the amount of H° inside the He + Stromgren sphere is negligible and in what follows 
we will adopt an ICF(Re) = 1.000. 

In Table 12 we present the helium abundance by mass Y(SMC), derived from region A. The 
F(SMC) values were derived from Table 9, the He ++ /H + value (that amounts to 2.2 x 10~ 4 for 
the T e range present in Table 12), and the Z values presented in Table 13. The Y(SMC) error 
presented in Table 13 is based only on the results of the self consistent method; by considering the 
T e (Bac) measurement the Y(SMC) error diminishes from 0.0018 to 0.0013 (and Y(SMC) diminishes 
by 0.0001). 

In Table 13 we present the helium, heavy elements, and oxygen abundance by mass of NGC 
346. The estimated error for O amounts to 0.06 dex. The Z value was obtained by assuming that 
O comprises 54% of the heavy elements by mass, this fraction was estimated by Carigi et al. (1995) 
from a group of ten irregular galaxies that includes the SMC. The estimated error for the Z value 
amounts to 0.08 dex. 



6. The Primordial Helium Abundance 

To determine the Y p value from the SMC it is necessary to estimate the fraction of helium 
present in the interstellar medium produced by galactic chemical evolution. We will assume that 

Y p = Y(SMC) - O(SMC)^. (26) 

In a recent review Peimbert & Peimbert (2000) derive Y v for the SMC from a similar discussion 
but using Z instead of O in equation (26); since the error in the O determination is smaller than 
in the Z determination (see Table 13) it is better to use O to determine Y p . 

To estimate AY/ AO we will consider three observational determinations and a few determi- 
nations predicted by chemical evolution models. 

Peimbert, Torres-Peimbert, & Ruiz (1992) and Esteban et al. (1999) found that Y = 0.2797 ± 
0.0060 and O = 0.0083 ± 0.0012 for the Galactic H n region M17, where we have added 0.08 dex 
to the oxygen gaseous abundance to take into account the fraction of these elements embedded in 
dust grains (Esteban et al. 1998). By comparing the Y and O values of M17 with those of region A 
we obtain AY/ AO = 5.45 ± 1.10. M17 is the best H n region to determine the helium abundance 
because among the brightest galactic H n regions it is the one with the highest degree of ionization 
and consequently with the smallest correction for the presence of He (i.e. 7CF(He) is very close 
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to unity e. g. Peimbert et al. 1992; Deharveng et al. 2000). It can be argued that the M17 AY/ AO 
value is not representative of irregular galaxies, because the yields are heavy element dependent 
and the O value is considerably higher than that of the SMC. 

From a group of 10 irregular and blue compact galaxies, that includes the LMC and the SMC, 
Carigi et al. (1995) found AY/ AO = 4.48 ± 1.02, where they added 0.2 dex to the O/H abundance 
ratios derived from the nebular data to take into account the temperature structure of the H n 
regions and the fraction of O embedded in dust; moreover they also estimated that O constitutes 
54% of the Z value. Izotov & Thuan (1998) from a group of 45 supergiant H n regions of low 
metalicity derived a AY/AZ = 2.3 ± 1.0; we find from their data that AY/AZ = 1.46 ± 0.60 by 
adding 0.2 dex to the O abundances to take into account the temperature structure of the H n 
regions and the fraction of O embedded in dust; furthermore from their data we also find that 
AY/ AO = 2.7 ± 1.2 by assuming that O constitutes 54% of the Z value. 

Based on their two-infall model for the chemical evolution of the Galaxy Chiappini, Matteucci, 
& Gratton (1997) find AY/ AO = 3.15 for the solar vicinity. Carigi (2000) computed chemical 
evolution models for the Galactic disk, under an inside-out formation scenario, based on different 
combinations of seven sets of stellar yields by different authors; the AY/ AO spread predicted by 
her models is in the 2.9 to 4.6 range for the Galactocentric distance of M17 (5.9 kpc). 

Carigi et al. (1995), based on yields by Maeder (1992), computed closed box models adequate 
for irregular galaxies, like the SMC, and obtained AY/ AO = 2.95. They also computed models 
with galactic outflows of well mixed material, that yielded AY/ AO values similar to those of the 
closed box models, and models with galactic outflows of O-rich material that yielded values higher 
than 2.95. The maximum AY /AO value that can be obtained with models of O-rich outflows, 
without entering into contradiction with the C/O and (Z— C — 0)/0 observational constraints, 
amounts to 3.5. 

Carigi, Colin, & Peimbert (1999), based on yields by Woosley, Langer, & Weaver (1993) and 
Woosley h Weaver (1995), computed chemical evolution models for irregular galaxies also, like 
the SMC, and found very similar values for closed box models with bursting star formation and 
constant star formation rates that amounted to AY / AO = 4.2. The models with O-rich outflows 
can increase AY /AO, but they predict higher C/O ratios than observed. 

From the previous discussion it follows that AY/ AO = 3.5 ± 0.9 is a representative value for 
models and observations of irregular galaxies. 

The Y p values in Table 14 were computed by adopting AY/ AO = 3.5 ± 0.9. The differences 
between Tables 12 and 14 depend on T e because the lower the T e value the higher the O value for 
the SMC. In Figure 5 we present our Y p value as well as the theoretical Y p value derived from Big 
Bang nucleosynthesis computations by Copi, Schramm, &: Turner (1995) for three neutrino species 
as a function of r], the baryon to photon ratio. 
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7. Conclusions 

The Y p value derived by us is significantly smaller than the value derived by Izotov & Thuan 

(1998) from the Y versus O/H linear regression for a sample of 45 BCGs, and by Izotov et al. (1999) 
from the average for the two most metal deficient galaxies known (I Zw 18 and SBS 0335-052), 
that amount to 0.2443 ± 0.0015 and 0.2452 ± 0.0015 respectively (see Figure 5). 

The difference could be due to systematic effects in the abundance determinations. There are 
two systematic effects not considered by Izotov and collaborators that we did take into account, 
the possible presence of H° inside the He + region and the use of a lower temperature than that 
provided by the [O m] lines. We consider the first effect to be a minor one and the second to be a 
mayor one but both should be estimated for each object. For further discussion of the first effect 
see the papers by Viegas et al. (2000) and Ballantyne et al. (2000). The second effect was first 
mentioned by Peimbert & Costero (1969, see also Peimbert 1995). 

From constant density chemically homogeneous models computed with CLOUDY (Ferland 
1996; Ferland et al. 1998) we estimate that the maximum temperature that should be used to 
determine the helium abundance should be 5% smaller than T e (0 III). Moreover, if in addition 
to photoionization there is additional energy injected to the H n region T e (He n) should be even 
smaller. 

Luridiana, Peimbert, & Leitherer (1999) produced a detailed photoionized model of NGC 2363. 
For the slit used by Izotov, Thuan, & Lipovetsky (1997) they find an JCF(He) = 0.993; moreover 
they also find that the T e (0 in) predicted by the model is considerably smaller than observed. From 
the data of Izotov et al. (1997) for NGC 2363, adopting a T e (He II) 10% smaller than T e (0 III) and 
AY/ AO = 3.5 ± 0.9 we find that Y p = 0.234 ± 0.006. 

Similarly, Stasinska & Schaerer (1999) produced a detailed model of I Zw 18 and find that the 
photoionized model predicts a T e (0 in) value 15% smaller than observed, on the other hand their 
model predicts an ICF(Re) = 1.00. From the observations of AA 5876 and 6678 by Izotov et al. 

(1999) of I Zw 18, and adopting a T e (He n) 10% smaller than T e (0 in) we obtain Y p = 0.237±0.007; 
for a T e (He II) 15% smaller than T e (0 III) we obtain Y p = 0.234 ± 0.007, both results in good 
agreement with our determination based on the SMC. Further discussion of these issues is presented 
elsewhere (Peimbert et al. 2000). 

The primordial helium abundance by mass of 0.2345 ± 0.0026(ler) — based on the SMC - 
combined with the Big Bang nucleosynthesis computations by Copi et al. (1995) for three light 
neutrino species implies that, at the one-u confidence level, Vi^h 2 is in the 0.0065 to 0.0091 range. 
For h = 0.65 the Y p value corresponds to 0.015 < fit < 0.022, a value considerably smaller than that 
derived from the pregalactic deuterium abundance, D p , determined by Buries & Tytler (1998) that 
corresponds to 0.041 < Clf, < 0.047(l<r) for h = 0.65. Our value is in very good agreement with 
the low redshift estimate of the global budget of baryons by Fukugita, Hogan, & Peebles (1998) 
who find 0.015 < $^5 < 0.030(ler) for h = 0.65 and is consistent with their minimum to maximum 
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range for redshift z = 3 that amounts to 0.012 < < 0.070. The discrepancy between Y p and D p 
needs to be studied further. 

In addition to the relevance for cosmology an accurate Y p value permits to determine AY/ AO, 
ratio that provides a strong constraint for the models of chemical evolution of galaxies (see for 
example Carigi et al. 1995, 1999). 

To increase the accuracy of the Y p determinations we need observations of very high quality 
of as many He 1 lines as possible to derive T e (He 11), iV e (He 11), T3889, and iV(He + )/iV(H + ) self- 
consistently. We also need observations with high spatial resolution to estimate the ICF(He) along 
different lines of sight. 

We are grateful to S. Torres-Peimbert for her participation on the observations presented in 
this paper. It is a pleasure to acknowledge several fruitful discussions on this subject with: L. 
Carigi, V. Luridiana, B. E. J. Pagel, M. Peha, E. Skillman, G. Steigman, S. Torres-Peimbert, and 
S. Viegas. We are also grateful to the referee for some excellent suggestions. MP received partial 
support from CONACyT grant 25451-E, MTR received partial support from a Catedra presidencial 
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Table 1. Positions, Sizes, and Observed Fluxes 



Region 


Position* 


Length b 


F(K8) 


Star c 


Spectral 




o 


,5 




-2 -1 
erg cm s 




Type 


1 


102.2 W 


0.0 


O A if! 

34.1 


6.23-14 






2 


58.0 W 


0.0 


22.4 


1.14-13 






3 


24.9 W 


0.0 


13.6 


1.63-13 






4 


2.1 E 


0.0 


18.6 


2.37-13 


355 


03V((f* )) 


5 


86.4 E 


0.0 


43.6 


2.54-13 






11 


123.9 W 


5.5 S 


24.8 


4.33-14 






12 


77.6 W 


5.5 S 


32.9 


9.48-14 






13 


42.5 W 


5.5 S 


25.8 


2.52-13 






14 


15.5 W 


5.5 S 


5.6 


4.51-14 


324 


04V((f)) 


15 


7.8 E 


5.5 S 


9.0 


1.08-13 


396 


07V 


16 


22.6 E 


5.5 S 


9.3 


1.39-13 


470+476 


08-08.5III+? 


17 


66.9 E 


5.5 S 


29.2 


2.22-13 






18 


92.4 E 


5.5 S 


13.1 


5.70-14 






A d 








11.6-13 






B c 








17.9-13 







a Relative to star 355: a = 00 h 57 m 19.82 s ; 5 = -72°26'40.4" (1950). 

b Slit is 1.6" wide and is oriented E-W. 

c Star numbers and spectral types from Massey et al. 1989. 

d Region A is the sum of Regions 2, 3, 5, 12, 13, 17, 18. 
e Region B is the sum of Regions 1 to 5 and 11 to 18. 
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Table 2. Line Intensities for Regions A and B 



A Id. Region A Region B a 



log F(\)/F(H (3) log I(X)/I(HP) 



3614 


He I 


-2.411 


-2.361 ± 0.040 




3634 


He I 


-2.475 


-2.425 ± 0.040 




3726+3729 


[0 ii] 


-0.049 


-0.003 ± 0.005 


+0.058 ± 0.005 


3750 


H12 


-1.546 


-1.501 ± 0.010 




3771 


Hll 


-1.421 


-1.376 ± 0.010 




3798 


H10 


-1.307 


-1.264 ± 0.010 




3835 


H9 


-1.176 


-1.134 ± 0.007 


-1.677 ± 0.015 


3o69 


fAT~ tttI 

[JNe inj 


— U.4U4 


— 0.OO4 ± 0.0U3 


— 0.3o6 ± 0.003 


3889+3889 


He I+H8 


-0.740 


-0.700 ± 0.004 


-0.915 ± 0.006 


3889 


He I 




-1.027 ± 0.008 




3967+3970 


[Ne m]+H7 


-0.568 


-0.533 ± 0.005 


-0.719 ± 0.006 


4026 


He I 


-1.767 


-1.733 ± 0.015 




4069+4076 


[S ii] 


-2.060 


-2.029 ± 0.020 




4102 


H5 


-0.621 


-0.591 ± 0.003 


-0.702 ± 0.005 


4340 


H 7 


-0.342 


-0.322 ± 0.003 


-0.400 ± 0.004 


4363 


[O in] 


-1.139 


-1.119 ± 0.005 


-1.175 ± 0.006 


4388 


He I 


-2.351 


-2.332 ± 0.020 




A A 71 

44 ( 1 


TT T 

lie I 


1 A QO 

— 1.4oz 


— 1.41b ± U.UUb 


— l.DO/ ± U.UUo 


4658 


[Fe in] 


-2.688 


-2.681 ± 0.050 




4686 


He ii 


-2.581 


-2.574 ± 0.040 




4711+4713 


[Ar iv]+He I 


-1.856 


-1.851 ± 0.015 


-2.028 ± 0.020 


4711 


[Ar iv] 




-2.029 ± 0.025 




4740 


[Ar iv] 


-2.242 


-2.238 ± 0.025 


-2.426 ± 0.030 


4861 


H/3 


0.000 


0.000 ± 0.000 


-0.023 ± 0.003 


4922 


He I 


-1.998 


-2.000 ± 0.008 


-2.310 ± 0.030 


4959 


[O in] 


+0.267 


+0.264 ± 0.002 


+0.238 ± 0.003 


5007 


[O in] 


+0.739 


+0.735 ± 0.002 


+0.707 ± 0.003 


5876 


He I 


-0.942 


-0.973 ± 0.005 


-1.022 ± 0.006 


6300 


[O i] 


< -2.016 b 


< -2.060 b 


< -2.083 b 


6312 


[S in] 


-1.727 


-1.771 ± 0.010 


-1.870 ± 0.012 


6563 


Ha 


+0.501 


+0.451 ± 0.002 


+0.451 ± 0.000' 


6584 


[N ii] 


-1.371 


-1.422 ± 0.005 


-1.365 ± 0.005 


6678 


He I 


-1.474 


-1.528 ± 0.003 


-1.582 ± 0.005 
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Table 2 — Continued 



A 


Id. 


Reg 


ion A 


Region B a 






log F{X)/F{H0) 


log I(X)/I(H (3) 




6716 


[S ii] 


-1.094 


-1.149 ± 0.003 


-1.142 ± 0.004 


6731 


[S ii] 


-1.237 


-1.293 ± 0.004 


-1.274 ± 0.004 


7065 


He I 


-1.615 


-1.675 ± 0.004 


-1.696 ± 0.005 


7136 


[Ar in] 


-1.065 


-1.127 ± 0.003 


-1.129 ± 0.003 


7281 


He I 


-2.135 


-2.199 ± 0.020 


-2.295 ± 0.025 


7320+7330 


[0 ii] 


-1.549 


-1.614 ± 0.005 


-1.511 ± 0.007 



a Given in log J(A)/J(Ha) + 0.451 (see text). 
b Upper limit. 

c Corrected for underlying absorption. 
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Table 3. Emission Line Intensities 21 
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u 




-L 
IE 


n 008 

u.uuo 


o 

u 


.319 


± 


o 

u. 


004 

UU^i 


o 

— u 


323 


± 


o 

u 


00^ 

UUO 


o 

— u. 


434 


-L 

ZLI 


004 

U.UU4: 


o 

u 


390 


-L 

zn 


o 

u 


004 

.UU4: 


^ouo 


1 

1 . 


1 ^ 


-1- 

ZL 


n 090 

u.uzu 


1 

1 


.152 


± 



u 


ni ^ 

UIO 


1 


.172 


± 


o 

u 


01 9 


1 

1 . 


1 71 


_j_ 


010 

U.U1U 


1 


039 


-L 

ZIZ 


o 

u 


007 

UU 1 




1 

1 . 


OO^i 


-L 

ZLI 


n o^o 

u.uou 


1 

1 


.404 


± 



u. 


090 


1 


,409 


± 


o 

u 


01 ^ 

UIO 


1 

1 . 


4RQ 
^iuy 


-L 

ZLI 


019 

U.U1Z 


1 

1 


499 


-L 


o 

u 


01 ^ 

-UIO 


4686 


< -2. 


.230' 


1 




< -2 


,700 d 






< -2. 


,873 d 






< -1. 


.974' 


1 




-2 


.199 


± 





.040 


4711+4713 










-2 


.275 


± 


0. 


.070 


-1 


,950 


± 


0. 


.035 










-1 


.697 


± 





.020 


4711 




















-2. 


188 


± 


0. 


.050 










-1 


.814 


± 





.030 


4740 




















-2. 


498 


± 


0. 


.080 










-2 


.035 


± 





.035 


4861 


0. 


.000 


± 


0.000 





.000 


± 


0. 


.000 





,000 


± 





.000 


-0. 


.045 


i 


0.003 





.000 


± 





.000 


4922 


-2. 


.019 


i 


0.035 


-2 


.004 


± 


0. 


.025 


-1. 


,968 


± 





.020 










-2 


.005 


± 





.015 


4959 


+0. 


.173 


i 


0.004 


+0 


.209 


± 


0. 


.003 


+0 


,264 


± 


0. 


.002 


+0. 


232 




0.003 


+0 


.311 


± 





.002 


5007 


+0 


.646 


i 


0.004 


+0 


.679 


± 


0. 


.003 


+0 


.736 


± 





.002 


+0. 


.704 


± 


0.003 


+0 


.783 


± 





.002 


5876 










-0 


.998 


± 


0. 


.020 


-0 


.962 


± 





.009 


-1. 


.029 




0.007 


-0 


.970 


± 





.007 


6312 


-1. 


.736 


i 


0.050 


-1 


.694 


± 





.030 


-1. 


.823 


± 


0. 


.025 


-1. 


.997 


± 


0.030 


-1 


.970 


± 





.030 


6563 


+0. 


.453 


± 


0.004 


+0 


.453 


± 





003 


+0. 


.453 


± 


0. 


.002 


+0. 


.451 


± 


0.000 c 


+0 


.447 


± 





.002 


6584 


-1. 


.227 




0.015 


-1 


.211 


± 


0. 


.010 


-1 


.705 


± 





.015 


-1. 


.485 




0.007 


-1 


.574 


± 





.008 


6678 


-1. 


.537 


i 


0.015 


-1 


.509 


± 


0. 


.012 


-1. 


.532 


± 


0. 


.009 


-1. 


.576 


i 


0.007 


-1 


.537 


± 





.006 


6716 


-0. 


.961 


± 


0.008 


-0 


.966 


± 


0. 


.006 


-1 


.526 


± 


0. 


.009 


-1. 


274 


± 


0.005 


-1 


.177 


± 





.005 


6731 


-1. 


.108 




0.010 


-1 


.110 


± 


0. 


.007 


-1. 


.667 


± 


0. 


.012 


-1. 


.430 




0.006 


-1 


.319 


± 





.005 


7065 


-1. 


.645 




0.020 


-1 


.643 


± 


0. 


.015 


-1. 


679 


± 





.012 


-1. 


.688 




0.009 


-1 


.673 


± 





.008 


7136 


-1. 


.152 


± 


0.010 


-1 


.090 


± 


0. 


.007 


-1. 


.115 


± 





.005 


-1. 


160 


± 


0.004 


-1 


.195 


± 





.005 


7320+7330 


-1. 


.410 




0.020 


-1 


.444 


± 


0. 


.015 


-1 


798 


± 


0. 


.015 


-1. 


.580 




0.010 


-1 


.873 


± 





.015 



a Given in log 7(A)/7(H/3). 
b Given in log J(A)/7(Ha) + 0.451 (see text). 
c Corrected for underlying absorption. 
d One sigma upper limit. 



Table 4. Emission Line Intensities' 



A 


11 




12 




13 








1 Ah 

14 




15 




16 




17 








18 






3726+3729 


+0.129 


± 0.010 


+0 


231 


± 0.007 


+0 


167 


± 


o 


006 


— 


115 +. 0.020 


+0 


051 ± 0.012 


+0 


328 zzz 0.008 


— () 


224 


± 


o 


008 


— 


430 


± 





020 


3835 


-1.174 


± 0.040 


_1 


155 


± 0.030 


— 1 


163 


+. 


o 


015 














— I 


155 


zzz 


() 


015 


— 1 


116 


± 





040 


3869 


-0.488 


± 0.012 


— 


453 


± 0.008 


— o 


422 


± 


o 


005 


— [) 


354 ± 0.020 


— 


494 ± 0.012 


— o 


485 zzz 0.012 


— () 


326 


± 


o 


005 


— o 


319 


± 


o 


012 


3889+3889 


-0.718 


± 0.015 


o 


707 


± 0.012 


_() 


704 


± 


o 


007 














— 


688 


± 


o 


010 


— 


698 


± 


o 


020 


3889 


-1.069 


± 0.030 


_ I 


043 


± 0.020 


_1 


037 


± 


o 


012 














_ I 


002 


± 


o 


012 


_1 


024 


± 


o 


035 


39674.397O 


-0.609 


± 0.015 


— 


555 


± 0.010 


— 


573 


± 


o 


006 














— () 


539 


± 


o 


006 


— 


537 


± 


o 


015 


4069+4076 






_ I 


930 


± 0.080 


_1 


988 


± 


o 


040 














_2 


131 


± 


o 


060 












4102 


-0.597 


± 0.015 


— 


593 


± 0.012 


— 


595 


± 


o 


006 






— 


982 ± 0.025 


— 


953 + 0.025 


— () 


589 


± 


o 


007 


— o 


583 


± 


o 


015 


4340 


-0.323 


± 0.008 


— 


332 


± 0.005 


— 


329 


± 


o 


003 


— 


509 ± 0.015 


— 


496 zL 0.009 


— 


405 + 007 

™o _L_ \J •\J\J 1 


— 


328 


± 


o 


003 


— 


325 


± 


o 


006 


4363 


-1.259 


± 0.025 


_1 


264 


± 0.015 


_1 


237 


± 


o 


010 


_1 


141 ± 0.040 


_1 


250 ± 0.025 


_1 


272 zL 0.025 


_1 


085 


± 


o 


010 


_ I 


022 


± 


o 


020 


4471 






_ 1 


410 


± 0.020 


_ 1 


419 


± 


o 


012 














_1 


415 


± 


o 


015 


_ 1 


428 


± 


o 


035 


4686 


< -2.181 d < 


2 


662 d 


^ 2 


815 d 


















^ 2 


433 d 






2 


609 d 






4711+4713 












-2 


152 


zzz 





035 














-1 


847 


± 





025 


-1 


694 


± 





050 


4711 












-2 


637 


zzz 





050 














-2 


023 


zzz 





035 


-1 


810 


± 





060 


4740 


































-2 


239 


zzz 





040 


-1 


965 


± 





070 


4861 


0.000 


± 0.000 





000 


± 0.000 





000 


± 





000 


-0 


077 ± 0.008 


-0 


071 ± 0.005 


-0 


025 ± 0.004 





000 


zzz 





000 





000 


± 





000 


4922 






-1 


995 


± 0.025 


-1 


982 


± 





012 














-2 


023 


zzz 





020 












4959 


+0.205 


± 0.004 


+0 


188 


± 0.003 


+0 


210 


± 





002 


+0 


282 ± 0.008 


+0 


151 ± 0.005 


+0 


135 ± 0.004 


+0 


294 


zzz 





002 


+0 


332 


± 





004 


5007 


+0.672 


± 0.004 


+0 


660 


± 0.003 


+0 


683 


± 





002 


+0 


757 ± 0.008 


+0 


628 ± 0.005 


+0 


597 ± 0.004 


+0 


767 


zzz 





002 


+0 


804 


± 





004 


5876 


-0.977 


± 0.020 


-0 


988 


± 0.012 


-0 


987 


± 





008 


-1 


074 ± 0.040 


-1 


047 ± 0.020 


-1 


205 ± 0.020 


-0 


983 


zzz 





008 


-0 


981 


± 





015 


6312 


-1.623 


± 0.040 


-1 


672 


± 0.035 


-1 


761 


± 





015 


-1 


953 ± 0.100 


-1 


931 ± 0.060 


-1 


961 ± 0.070 


-1 


869 


zzz 





025 


-1 


980 


± 





040 


6563 


+0.445 


± 0.004 


+0 


446 


± 0.003 


+0 


446 


i 





002 


+0 


451 ± 0.000° 


+0 


451 ± 0.000° 


+0 


451 ± 0.000° 


+0 


452 


zzz 





002 


+0 


443 


± 





004 


6584 


-1.381 


± 0.020 


-1 


228 


± 0.012 


-1 


273 


zzz 





008 


-1 


546 ± 0.040 


-1 


379 ± 0.020 


-1 


137 ± 0.012 


-1 


732 


zzz 





015 


-1 


769 


± 





025 


6678 


-1.546 


± 0.020 


-1 


542 


± 0.012 


-1 


532 


zzz 





007 


-1 


672 ± 0.040 


-1 


600 ± 0.020 


-1 


719 ± 0.020 


-1 


526 


zzz 





008 


-1 


546 


± 





020 


6716 


-1.118 


± 0.012 


-0 


966 


± 0.006 


-1 


041 


zzz 





004 


-1 


396 ± 0.030 


-1 


140 ± 0.012 


-0 


985 ± 0.008 


-1 


246 


zzz 





005 


-1 


512 


± 





020 


6731 


-1.306 


± 0.015 


-1 


124 


± 0.007 


-1 


184 


zzz 





005 


-1 


548 ± 0.035 


-1 


251 ± 0.015 


-1 


075 ± 0.009 


-1 


389 


zzz 





006 


-1 


677 


± 





020 


7065 


-1.627 


± 0.025 


-1 


690 


± 0.015 


-1 


687 


zzz 





010 


-1 


815 ± 0.050 


-1 


759 ± 0.025 


-1 


907 ± 0.030 


-1 


669 


zzz 





010 


-1 


679 


± 





020 


7136 


-1.097 


± 0.012 


-1 


098 


± 0.007 


-1 


096 


± 





004 


-1 


168 ± 0.020 


-1 


174 ± 0.012 


-1 


175 ± 0.010 


-1 


145 


± 





005 


-1 


202 


± 





012 


7320+7330 


-1.585 


± 0.025 


-1 


385 


± 0.015 


-1 


406 


zzz 





012 


-1 


645 ± 0.040 


-1 


486 ± 0.020 


-1 


189 ± 0.020 


-1 


864 


zzz 





015 


-1 


945 


± 





030 


a Givcn in 
k Given in 


log J(A)/J(H/3). 

log I(\)/I(Ha) + 0.451 (sco text) 















































c Corrected for underlying absorption. 



One sigma upper limit. 
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Table 5. Emission Equivalent Widths' 



Region 


Ha 


H/3 


H 7 


H<5 


H9 


4471 


4922 


5876 


6678 


7065 


i 
± 




oo 


9^ 


1 1 


9 9 


9 ^ 


n q 


10 


7 n 


U.(J 


2 


720 


190 


75 


38 


10 


5.5 


2.0 


28 


9.5 


8.5 


3 


1650 


380 


140 


65 


17 


12 


4.5 


60 


20 


16 


4 


115 


19 


5.0 


1.5 


-0.65 b 


0.55 


0.35 


2.8 


1.2 


1.1 


5 


785 


210 


90 


45 


11 


7.5 


2.2 


28 


8.5 


7.0 


11 


580 


140 


60 


30 


7.5 


5.0 




22 


7.0 


5.5 


12 


790 


200 


75 


38 


9.5 


7.2 


1.4 


33 


8.5 


7.0 


13 


1760 


350 


115 


60 


13 


11 


4.0 


65 


21 


16 


14 


55 


7.0 


2.0 










1.4 


0.50 


0.45 


15 


200 


22 


4.8 










4.5 


2.0 


1.9 


16 


69 


8.0 


2.5 






-0.55 b 




1.0 


0.29 


0.16 


17 


1950 


340 


120 


55 


14 


10 


3.5 


70 


21 


19 


18 


1800 


300 


100 


60 


13 


9.5 


4.0 


65 


20 


18 


A 


1250 


250 


99 


44 


12.5 


8.5 


2.5 


46 


14.7 


11.2 


B 


230 


35 


11.1 


4.7 


0.43 


0.7 


0.18 


6.5 


2.7 


2.2 



a In Angstroms. 
h EW in absorption. 
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Table 6. Temperatures and Densities 



Region T(0 III) T(0 n) N e {S n) 



1 


13 


730 


± 


250 


11 


650 


± 


450 


40 ± 


35 


2 


13 


350 


± 


190 


11 


150 


± 


400 


50 ± 


25 


3 


12 


500 


± 


130 


12 


600 


± 


350 


55 ± 


45 


4 


12 


840 


± 


130 


13 


600 


± 


350 


15 ± 


20 


5 


13 


540 


± 


100 


10 


950 


± 


250 


55 ± 


20 


11 


12 


230 


± 


250 


10 


500 


± 


600 


14 a ± 


50 


12 


12 


320 


± 


170 


12 


100 


± 


350 


15 ± 


25 


13 


12 


360 


± 


130 


12 


600 


± 


300 


50 ± 


20 


14 


12 


580 


± 


390 


13 


550 


± 


1150 


25 ± 


150 


15 


12 


800 


± 


270 


12 


550 


± 


550 


150 ± 


60 


16 


12 


890 


± 


310 


12 


550 


± 


600 


220 ± 


45 


17 


13 


120 


± 


140 


11 


100 


± 


300 


55 ± 


20 


18 


13 


440 


± 


280 


13 


500 


± 


950 


14 a ± 


65 


A 


13 


070 


± 


50 


11 


810 


± 


160 


50 ± 


15 


B 


12 


755 


± 


90 


12 


360 


± 


220 


85 ± 


20 



a Low density limit. 
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Table 7. Temperatures and Densities for Region A 



T e (°K) 


e 




7V e (cm- 3 ) 


(0 m)=13070 ± 100 






(rms)= 14 ± 3 


(0 n) = 11810 ± 300 






(S n) = 50 ± 15 


(0 n+0 m)=12860 ± 100 








(Bac)=11800±500 


(Bac, n+0 in) 


=0.023 ±0.011 




(He n) Bac =11890±500 








(He n) sc =11950±560 


(He II, n+0 m) 


=0.021 ±0.012 


(He n)sc=143±50 


((He n) } = 11920 ± 370 


(Bac, He n, O n+O m) 


=0.022 ±0.008 


((He n))=146±50 



Table 8. Ionic Abundances 51 for t 2 = 0.0013 b 



Region 


0+ 


0++ 


N+ 


Ne++ 


S+ 


S++ 


Ar++ 


Ar 3+ 


i 
i 


7 AA 


7 7Q 


^ 87 


u.yo 


^ zL9 


fill 






2 


7.47 


7.86 


5.91 


7.04 


5.44 


6.19 


5.61 




3 


7.10 


7.99 


5.45 


7.15 


4.93 


6.16 


5.64 


4.65 


4 


7.27 


7.93 


5.67 


7.12 


5.15 


5.95 


5.57 




5 


7.01 


7.94 


5.53 


7.11 


5.21 


5.89 


5.49 


4.98 


11 


7.51 


7.96 


5.81 


7.10 


5.35 


6.40 


5.68 




12 


7.60 


7.94 


5.96 


7.12 


5.50 


6.34 


5.67 




13 


7.53 


7.95 


5.91 


7.15 


5.42 


6.25 


5.66 


4.56 


14 


7.22 


8.01 


5.62 


7.19 


5.05 


6.02 


5.58 




15 


7.38 


7.86 


5.77 


7.03 


5.32 


6.02 


5.56 




16 


7.65 


7.82 


6.01 


7.03 


5.49 


5.98 


5.55 




17 


7.06 


7.97 


5.40 


7.16 


5.17 


6.05 


5.56 


4.81 


18 


6.83 


7.97 


5.34 


7.13 


4.88 


5.90 


5.49 


5.07 


A 


7.28 


7.94 


5.69 


7.13 


5.27 


6.15 


5.59 


4.81 


B 


7.41 


7.97 


5.81 


7.17 


5.33 


6.11 


5.63 


4.64 



a Given by log N(X)/N(H+) + 12. 

b We have adopted a two temperature model with T(0 II) and T(0 in), 
inside the + and ++ zones the temperature is uniform (t 2 = 0.000), 
but the t 2 over the entire model turns out to be 0.0013. Notice that our 
preferred t 2 value is 0.022, see Tables 7 and 11. 
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Table 9. 


, iV(He 4 


)/N(K 


f ) a and x 2 


for Re£ 


^ion A 


Te(K) 


53 


100 


N e (cm- 3 ) 
143 


162 


247 


11200 


805 
(83.2) 


798 
(47.7) 


793 
(26.4) 


791 
(20.0) 


781 
(8.24) b 


11800 


806 
(38.6) 


799 
(15.9) 


793 
(7.37) 


790 
(6.59) b 


780 
(20.4) 


11950 


806 
(30.8) 


799 
(11.7) 


793 

(6.53) b - c 


790 
(7.25) 


779 
(27.7) 


12400 


807 
(15.0) 


799 793 
(7.17) b (12.5) 


790 
(17.9) 


778 
(58.6) 


13000 


809 
(9.72) b 


800 
(18.2) 


793 
(38.4) 


790 
(50.2) 


777 
(118) 


a Given 


in units of 10 4 , 


X 2 values 


in parenthesis. 



b The minimum x 2 value at a given temperature. 

c The smallest x 2 value for all temperatures and den- 
sities, thus defining T e (He n) and N e (Re n). 
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Table 10. Ionic Helium Abundances a 



Region 


3889 


4471 


4922 


5876 


6678 


7065 


(He+/H+) 


c 


1 




838 


739 




778 


849 


802 


2.07 


2 


802 


800 


764 


755 


830 


853 


815 


2.21 


3 


828 


791 


831 


820 


787 


784 


804 


2.08 


5 


789 


768 


763 


806 


778 


796 


788 


0.525 


11 


735 






793 


762 


885 


794 


4.03 


12 


780 


789 


780 


763 


769 


765 


773 


3.10 


13 


791 


773 


804 


775 


787 


771 


783 


2.36 


17 


850 


783 


735 


786 


805 


788 


797 


0.873 


18 


815 


757 




785 


763 


785 


780 


0.719 


A 


809 


779 


772 


800 


795 


792 


793 b 


1.57 



a Given by 10 4 x N(Re + )/N(R + ). 
b Average of 9 lines, see Table 9. 
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Table 11. Total Abundances a 



Element 


SMC b 


NGC 346 c 
t 2 = 0.0013 f t 2 = 0.0220 s 


Sun d 


M17 e 


N 


6.46 


6.44 


6.51 


7.92 


7.90 





8.02 


8.07 


8.15 


8.83 


8.87 


Ne 


7.22 


7.22 


7.30 


8.08 


8.02 


S 


6.49 


6.50 


6.59 


7.33 


7.31 


Ar 


5.78 


5.74 


5.82 


6.40 


6.60 



a Given by logiV(X)/iV(H) + 12. 

b Dufour 1984. 
c This paper. 

d Grevesse & Sauval 1998. 

e Peimbert et al. 1992; Esteban et al. 1999. 

f Minimum t 2 , see first paragraph in section 4. 
g Preferred t 2 , see section 3 and Table 7. 
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Table 12. Y(SMC) 



T e (K) 






iV e (cm _3 ) 








53 


100 


143 


162 


247 


11200 


0.2431 


0.2416 


0.2404 


0.2399 


0.2377 a 


11800 


0.2435 


0.2419 


0.2405 


0.2399 a 


0.2375 


11950 


0.2436 


0.2420 


0.2405 a 


0.2399 


0.2374 


12400 


0.2439 


0.2421 £ 


1 0.2406 


0.2399 


0.2372 


13000 


0.2443 a 


0.2423 


0.2407 


0.2400 


0.2370 



a Entries that correspond to minimum \ 2 values, see 
Table 9. 
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Table 13. Chemical Composition* 1 



Element 




NGC 346 b 


M17 b ' c 




t 2 = 0.0013 t 2 = 0.0220 




Y 


0.2445 


0.24050±0.00180 


0.2797±0.0060 


Z 


0.00263 


0.00315±0.00063 


0.0212±0.0030 


o 


0.00142 


0.00171±0.00025 


0.0083±0.0012 


AY/AZ 




1.9±0.5 


2.13±0.5 


AY/AO 




3.5±0.9 


5.45±1.1 



a Given by mass. 
b This paper. 

c Peimbert et al. 1992; Esteban et al. 1999. 
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Table 14. Y p Derived from the SMC 



Te(K) 






iV e (cm _3 ) 








53 


100 


143 


162 


247 


11200 


0.2363 


0.2348 


0.2336 


0.2331 


0.2309 a 


11800 


0.2373 


0.2357 


0.2343 


0.2337 a 


0.2313 


11950 


0.2376 


0.2360 


0.2345 a 


0.2339 


0.2314 


12400 


0.2384 


0.2366 £ 


1 0.2351 


0.2344 


0.2317 


13000 


0.2395 a 


0.2375 


0.2359 


0.2352 


0.2322 



a Entries that correspond to minimum \ 2 values, see 
Table 9. 
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Fig. 1. — Blue spectrum for region A. 
Fig. 2. — Red spectrum for region A. 

Fig. 3. — Spectrum of region A near Ha that shows [N II], [S II], and He I lines. 

Fig. 4. — Spectra of NGC 346 with and without underlying absorption. The vertical scale is for 
the lower spectrum (region A). The flux of the upper spectrum (region B) was normalized to the 
Ha emmision line flux for the lower spectrum. 

Fig. 5. — Y p versus rj diagram from the Big Bang nucleosynthesis computations by Copi et al. 
(1995). The low Y p box is by us while the high Y p box is by Izotov et al. (1999). The dashed curves 
correspond to la error bars, the boxes include la errors in the observations and the theoretical 
cmputations added in quadrature. 
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